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Abstract 

A derivation of the neutrinoless double beta decay rate, specially adapted for the 
nuclear structure calculations, is presented. It is shown that the Fourier-Bessel expansion 
of the hadronic currents, jointly with the angular momentum recoupling, leads to very 
simple final expressions for the nuclear form factors. This greatly facilitates the theoretical 
estimate of the half life. Our approach does not require the closure approximation, which 
however can be implemented if desired. The method is exemplified for the BB decay 
48 Ca — > 48 Ti, both within the QRPA and a shell-model like model. 
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1 Introduction 



The standard model (SM) of electroweak interactions has been brilliantly confirmed 
by a host of experiments. But, about the properties of neutrinos it says very little. 
Actually, in the SM it is postulated that: i) the neutrinos are the only fermions without the 
right-handed partners, and ii) their masslessness is dictated by the global lepton-number 
symmetry, and not by a fundamental underlying principle, such as gauge invariance for 
the photon. More, whether neutrinos behave so "trivially" as required by the SM is one 
of the most fundamental open questions of the present-day physics. 

It has been known for a long time Jl|, |3], |], || that neutrinoless double beta decay 
{PPou) is a very sensitive probe of lepton number violating terms in the Lagrangian such as 
the Majorana mass of the light neutrinos, right-handed weak couplings as well as the Higgs 
exchange ||, right-handed weak coupling involving heavy Majorana neutrinos J7|, massless 
Majoron emission 0, |8|, |[| [H| [11), and R-parity breaking in the supersymmetric model 
fT2| , |l3|j . Thus, if the /3(3q u decay is someday observed experimentally it would hint new 
physics beyond the SM. But, even if it is not observed, the measured limits on its transition 



probability, which are steadily improving [TJj], could be translated into more stringent 
constraints on the parameters of the just mentioned new theoretical developments. 

Yet, the extraction of these constraints from the data is only possible when we know 
how to deal with the nuclear structure involved in the f3(3ou decay. This is not at all an 
easy task, because of: 

(1) the nuclear hamiltonian is only roughly known to the extent that the choice of the 
appropriate parametrization is an art, 

(2) there is in general a very large amount of nuclear states involved in the calculation, 

and 

(3) the formulas for the fifiov decay rate are rather complex and difficult to implement 
in a nuclear structure calculation. 

In the present work we derive simple expressions for the nuclear matrix elements, 
especially tailored for the nuclear structure calculations. The simplification mainly comes 
from the Fourier-Bessel expansion of the term exp[ik- (rj— r 2 )] in the transition amplitude, 
and in doing the integrations in the following order: first on dQ^, then on dr\ and dr2, 
and finally on k 2 dk [15]. So far, the same procedure has been applied for the evaluation 
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of the matrix elements M F , M GT |T7| and M R JTI| that arise from the electron s-wave. 



Here we also dealt with the p-wave matrix elements that are relevant when the admixture 
of the right-hand lepton current is considered. Other studies on the subject are those 
of Vergados et al. [|18| , who derived the formulas directly in the momentum space, and 
those of Suhonen, Khadkikar and Faessler JT^|, who worked in a framework of a relativistic 



quark confinement model. 

This paper is organized as follows: In sec. 2 we discuss the basic mechanism for the 
(3(3q v decay, presenting the effective hamiltonian and the transition amplitude in the form 
convenient for the multipole expansion, which is carried out in sec. 3. In sec. 4 we give 
the detailed formulas for the nuclear matrix elements and discuss the nuclear structure 
calculations involved in the problem. Summarizing conclusions are drawn in sec. 5. 

2 Effective Hamiltonian and the half life 

The 0u(3(3 half life 

[To„(0 + - 0+)]- 1 = g|, (1) 

for the decay from the state | /) in the (N, Z) nucleus to the state \f) in the (N — 2,Z + 2) 
nucleus (with energies Ej and Ep and spins and parities J n = + ), is evaluated via the 



second order Fermi's golden rule. Thus the decay rate (in % = c = m e units) is pO 



with 



T 0u = 2nY: I |ifr»(ei, 0.) ?S{e x + e 2 + E F - E t ) f 1 ^ , (2) 

V J ( 2n ) Ej- E N -ei-u 

where e = (e, p, s e ) [y = (u, k, s u )) stands for the energy, momentum and spin projection 
of the electron (neutrino), and n runs over all levels in the (N — 1,Z + 1) nucleus. 
The effective weak hamiltonian reads p], £|. || 

/"f In „ 

H * = ^ £/ dx[7Lft,(x) J£j(x) + j^(x)J^(x) + h.c), (4) 
where the summation goes on the number on lepton generations, 

f LM {*) = 2^(x)7^N L ^(x); P L>/? = I(1 T75 ), (5) 
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are the leptonic currents, formed out from the electron field ^(x) and the Majorana 
neutrino field N^(x) of mass mi, and 

J£{(x) = V*J* (x), Jg(x) = K,(AJ^( X ) + (x)), (6) 

contain the hadronic (V^fA) currents J^r- ^ U e e and V e £ are the neutrino mixing matrices 
for the left- and right-handed sectors, and A and rj are the strengths of admixtures of the 
(V + A) current. 

Within the non-relativistic impulse approximation the hadronic currents read, 



JlA*) = ( v Pv(x)TP. 4 (x),jv(x) T j. 4 (x)j (7) 

where @ 

Pv(x) = # v ^r+5(x- r„), 

n 

Pa(x) = T^TrH^i^n ■p„5(x-r n ) + 5(x-r n )<T n -p n ], (8) 



2M N n 



9a 
9v_ 

j A (x) = g A ^T+a n 5(x-r n ), 



jv(x) = — 7 H r n + W( x - r J + - r„)p„ + /» Vxtr^(x - r, 

Z1V1 N n 



are the one-body vector (V) and axial-vector (A) densities and currents, M N is nucleon 
mass and f w = 4.7 is the effective weak-magnetism coupling constant. 
Merging (Q) into (0) and performing the s v -summation one gets [|[: 

C 2 2n r r rlh 

Rou = ^=EE E / ^rfy / ™(F|J,1(y)e^|v)(v|^](x)e-*- x |/) 

l=l N a ,(3=L,R J J y Z71 )' 

v r -i P / m jfa, y>T^ - k ■ 7 + m £ ) P a7 ^ c (e 1 , x) 

x l — ^ei,e 2 J 7 ; ttt , l y J 

w(ei + u + h N - hjj) 

where ip(ei,x) and ■0(e2,x) are the wave functions of the emitted electrons, and the 
operator P (61,62) interchanges the particles 6\ and 62- The structure of the eq. (0) 



1 We do not consider the admixture of the hadronic (V + A) current into j£, since its contribution to 
the (3/3 decay amplitudes is negligible Q. 

2 See eq. (3D-18) in ref. pjj . The correspondence between the non-relativistic approximations used 
here and that prevailingly employed in the studies of the pPov decay [| ||, 0, can be find on p. 516 
of the Walecka's book p|. 
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suggests that it might be convenient to introduce the Fourier transforms of the quantities 
defined in (§), i-e., 

p(k) = |rfxp(x)e- kx , 

j(k) = /rfxj(x)e- jk - x . (10) 



Next, ensuing the usual procedure |2|, f|, ||, we evaluate the sy2 and pi/ 2 contributions 
of the electron wave functions to the amplitude Rq v . The first ones give rise to the 
following k and n dependent nuclear moments 

M F (k,w) = (F\ Pv (-k)\N}(N\p v (k)\l), 
M GT (k,N) = (F|j A (-k)|AT).(iv|j A (k)|z), (11) 

M R (k, N ) = -iRk- (F|j A (-k)|iv)x(jv|j v (k)|j), 
where R is the nuclear radius, and the second one to 
M F ,(k,N) = 2v / 3^F|p(°)(-k)|^)(7v|p v (k)|/), 

M GT ,(k,7v) = 6i( F \jW(-k)\N)-{ N \j A (k)\i), (12) 
M T (k,N) = 2i(Fbl 21) (-k)|iv)-(Ar|j A (k)|/}, 

M P (k,N) = -V2i[(F\jW(-k)\N)(N\ Pv (k)\l) - (F\U-k)\N) ■ Hp«(k)|/)], 

where we have introduced the tensor operators 

p (J) (k) = |rfxp(x)(k®x)( J ) e - ik - x , 



3 (LJ) 



(k) = LJ- 1 /dx[j(x)® (k®x)W]^e- ik - x , (13) 



with L = \J2L + 1. The explicit form of the matrix elements defined in ( [TT]) and ( |T2"D are 
shown in the appendix A. 

We now define the nuclear matrix elements 



Aug 



J dkv (k,u N )M x (k, n) for x = f, gt, f',gt', p, _r,t, (14) 
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and 

R 



dkv u (k,uj N )M x (k,N) for x = f 1 gt ) (15) 



^9 A N 

with 

^(fc,^) = ^(fc,^ jV ) = - (16) 

and 

= ^ - \ (E, + E F ) . (17) 

In deriving the expression ( |T7| ) we have approximated the electron energies as 6^2 — 
(£J Z — E F )/2. We have also neglected the neutrino mass in comparison with k, i.e., we 
have taken uo = k. 

For the transition amplitude we get 

g 2 G 2 5 

Rov{ei, e 2 ) = A ^ R ^ J2 z kL k (ei, e 2 ), (18) 



where 



Z x = <m v > (M F - Mar), 

Z 2 = <r]> (Mgtuj + M Fu) )+ < A > (M FuJ - Mgtuj), 

Z 3 = A< V >M R , (19) 

Z 4 = \i[< A > (M' GT - 6M T + ZM'p)— < r) > (M' GT - 6M T - 3M' F )], 

Z 5 = Ai<r]> M P , 

encompass the hadronic matrix elements, as well as the parameters 

< m v > = J2' m t U ee, 

t 

<X> = Xj2'U e eV ee , (20) 

i 

<T]> = rf^UeiVet, 

i 

where the summation J2/ goes only on the light neutrinos j|, [5]. The leptonic matrix 
elements L\ (61,62) are displayed in the appendix B. 
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Finally, by performing the integrations (summations) on the electron states indicated 
in (Q), we obtain the familiar expression for the half life [0, [|] 



[T 0u (0 + -> 0+)]- 1 = <m„ > 2 C,+ < \ > 2 C 2 + < V > 2 C 3 (21) 
+ < >< A > C 4 + < m v >< r\ > C 5 + < A >< rj > C 6 , 

where 

d = (M F -M GT ) 2 £ 1; 

c 2 = MtG2 + ^M 2 + g 4 -^M 2 _M 1+ g 3 , 

c 3 = Mi + g 2 + -Ml_Qt - -m 2+ M!_^3 + M 2 R g 9 + M R M P g 7 + M 2 g 8 , 

C A = (M F -M GT )[M 2 _g 3 -M 1+ g 4 ] } (22) 

c 5 = -(M F -M GT )[M 2+ g 3 -M 1 _g i + M R g fi + M P g 5 ], 

C 6 = -2M 2 _M 2+ g 2 + ~ [M 2 „Mi_ + M 2+ M 1+ ] £ 3 - ~Mi_Mi+0 4 , 
contain the usual combinations of the matrix elements 

M 1± = M' GT - 6M T ± 3M' F , 

M 2± = M GTu1 ±M Fuj -\m 1t , (23) 

y 

and the kinematical factors 

g A A G 4 ( 2tt«Z ^ 2 



32R% 5 ln 2 VI- e 



-1-notZ 



MTo). (24) 



The electron phase-space factors jF fc (T ), as a function of the maximum kinetic energy 
T = Ej — E F — 2, are listed in the appendix B. 

It might be important to stress that, within the procedure followed here to derive the 
result fll8|), we do not need to recur at all to so called closure approximation (CA). (Remind 
that the CA implies: i) to supplant the energies E N by an average values < E N >, and 
ii) to use the closure relation J2 N \ N )( N \ = 1 f°r the intermediate states.) When reworked 
in the CA, the moments ( |i~4"|) and ( |i~5|) are directly comparable with those that appear in 
the literature § g § (§ §. 
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3 Multipole expansion and angular momentum 

recoupling 



The starting point for the multipole expansion of the hadronic current is to use the 
Fourier-Bessel relation 

e* k ' r = A^ L JL{kr){Y L {^)-Y L {v)) 

L 

= 47r^2 L (-l) L Lj i (A;r)[r L (k)®F L (f)] , (25) 

L 

in the equations exhibited in the appendix A. Then we perform the angular momentum 
recoupling, and rewrite the nuclear moments (|ll|) and fll2|) in terms of the one-body 
spherical tensor operators 

YajmO) = ^2r+r^ x (kr n )Y JM (r n ), 

n 

SlLJAlW = Y, T n r nJx( kr n)[(Tn®Y L (r rl )}j M , 

n 

P LJ M{k) = Y, T nJL{kr n )[p n ®Y L {r n )] JM . (26) 

n 

Finally, the angular integration on dQ^ is done. We illustrate the procedure by sketching 
in the appendix C, the derivation of a part of the final formula for the nuclear matrix 
element M R . Proceeding in a similar way with the remaining matrix elements we obtain: 

M F = 4vrR^] ^Jv(k,uj N )k 2 dk(F\YUk)\N)-(N\YUk)\i), (27) 



Jn 



M GT = 4vrR^(-l) 1+L+J |t;(A ; ,c, JV )A; 2 ^(F|S LL j(A:)K)-^|S LL j(A:)|/), (2f 



LJn 

9v \ ST^ ;L-J+1 



M F , = -8vrR — J2 lL ~ J+1 ( J1 \L)(Jl\L) 



A 



LJn 



x /i;(fc,a; JV )fc 3 dfc(F|Yi t7 (fc)| J v)-( J v|Y° J (fc)|/), (29) 



M GT > = 8ttR Y, i L ~ L ' +1 (-l) L ' +J (L'l\L)(L'l\L) 

LL'Jn 

X / v(k,U N )k 3 dk(F\S{ LIJ (k)\N) ■ (N\S L>UJ (k)\ F ), (30) 
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U« = E i L+L \~l) J fv(k,u; N )^dk( F \Sl LJ (k)\N) 

M N 9a LUJn ■' 



(N\S° L , L ,j(k)\l) 



2V6(-1) L+J L{^ \ \, \ (Ll\L>)(N\P L ,j(k)\i)\ , (31) 



8 LL ,-{J1\L){J1\L') 

L J 1 
1 1 V 



M T = 40ttR i L+L ' +1 L 2 (lL\J')(lL\L') 

LL'J'Jn 

x f v{k^ N )k*dk{F\S l LJIJ {k)\N) ■ (N\S° L , LfJ (k)\i), (32) 





2 






2 






L 






J 





M P = 8irVQR— T i L+J+l j(Jl\L)(Jl\L){ f f * 1 

x [vik^^dkiFlSljMl^-iNlYUkM, (33) 



where (J1|J) is a short notation for the Clebsh- Gordon coefficient (L010| JO). The for- 
mulas for the matrix elements M FLU and M GTu) are obtained from those for M F and M GT 
with the replacement v(k,u N ) — ► v w {k,uj N ) (see eqs. ( |14D and (|15|)). 
The evaluation of the Qv(3(3 matrix elements encompasses: 

i) the appraisal of the scalar product 

(F\Jj(k)\N)-(N\Tj(k)\l), (34) 

where Tj(k) represents any of the one-body operators displayed in fl26D , and 

ii) the integration on the neutrino momentum k. 

More details on these two steps are given in the next section. 
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4 Nuclear structure calculations 



To evaluate the matrix elements (B3j) it is convenient to rewrite the operators (|26| ) 
from the Hilbert space to the Fock space ET[, i.e., 



Tj M (k) = J~ ll £(p\\Tj(k)\\n) (a\a n ) 

pn 



JM 



(35) 



In this way we get 

-£(F\Tj(k)\N) ■ (N\Tj(k)\i) = ^ (0-}\Tj(k)\JZM) ■ (J:M\Tj(k)\0f) 



M 



carM 



E (PllTAm^ff^pnp'n'^^ip'llTjmn'}., 



where 



ff h (pnp'n';JZ) = .Z" 2 (0+ 1 1 « ^) J7r 1 1 1 1 (a^a^ ) ^ 1 1 0+) , 



(36) 



(37) 



J f I I \ p ft / j " i i "a / \" a 1 1 \—p 

is a two-body state dependent particle-hole (ph) density matrix, and the index a labels 
different intermediate states with the same spin J and parity tt. 

Within the CA we can sum over a, and deal with the state independent ph density 
matrix 



fg(pnp'ri; r) = XyWn'; J£) = ^(O+l [(aJo fi )j*(a^a fi Oj*] |0 f +). 
which is related with the particle-particle (pp) density matrix 

ffP(pp'nn'; J*) = J'^Q+l [(ajaj,)^(o^)^] |0+), 
by a Pandya like relation 

° T \ P^(pp'nn';r). 

o> 1 J 

The reduced single-particle pn form factors for the one-body operators defined in 



pff (pnp'n'; J*) = £(-)i*+V+^ { ^ J« J 



(39) 



(40) 



arc 



(pllYJjWIln) 
(p\\Sl LJ (k)\\n) 
(p\\P L j(k)\\n) 



(4tt) 2^ J0J (pn)^(pn;A;), 
(47r)-^ L1J (pn)^(pn;fc), 



(41) 



>(-) 



(+)/ 



?(+)< 
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with the angular parts! 

W LS j(pn) 
Wif(pn) - 

and the radial parts 

Rl{pn;k) = 



lp 2 fa 



\/2SJLi n j n j p (l n L\l p ) I L S J 

T i(-l)^+J+hjti p j p j n (ln + \ T \)^{lpL\ln T 1) 

| lp h \\\ L J 1 } , 
1 jn In J I 1 In ^n ~F 1 4 



Rl(l p ,n p ,l n ,n n ;k) = I u np ^{r)u n ^i n {r)j L {kr)r 2+K dr 



(42) 



' d 2Z n + 1 ± r 
dr 2r 



u nn ,i n {r)j L {qr)r 2 dr. 



(43) 



To carry out the numerical calculation of eqs. (27) - (33) it is convenient to group 
separately the angular and the radial parts. For instance, M GT can be cast in the form 

Mgt = -E(") L E ^(pnpV^^MjH^ijfpVl^fwV;^), (44) 

LJ™ pp'nn' 

where the two-body radial integrals are defined as 

ni L ,(pnp'n'; u JS ) = rJ dkk 2+K v(k; u JS )R L (pn; k)R° L ,(p'n; k). (45) 

One manner to include the effects of the finite nucleon size (FNS) and the two-nucleon 
short-range correlations (SRC) on the 13(5q v moments has been explicated in ref. [IB], with 
the result: 

,2 \4 i 



v (k, 



A 2 + k 2 



Trkk, 



In 



+ 



2rckk r 



3 1 



,n ri 



c Ln=l 



xl + In 



X- 

Xj 



(46) 



where A = 850 MeV is the cutoff for the dipole form factor in the FNS correlations, 

A 2 

x± 



A 2 + (k± k c 



(47) 



3 We use here the angular momentum coupling 
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and k c = 3.93 fm -1 is roughly the Compton wavelength of the cu-meson in the SRC 
correlations. 

The integration on the neutrino momentum k is simplified when the harmonic oscillator 
radial wave functions are employed. Then the following relations among the one-body 
radial integrals are valid: 

Rl{pn;k) = (2i/)-£ {(2l n + 2n n + 3)?R° L (k;l p ,n p ,l n + l,n n ) 

- (2n n )^R° L (k;lp,n p ,l n + l,n n -l)j, 
i 

Rf\pn]k) = ± (j0 2 {(2Z n + 2n n + 2=Fl)3 J R°(fc;/ p ,n p ,Z n =Fl,n n ) 

+ {2n n + l±l)$RP L (k;l p ,n p ,l n =Fl,n n ±l)}, (48) 

where v = Muj/Ti is the oscillator parameter, and the ^-integration in the matrix elements 
( p7|) - ( j33[) only involves the radial integrals (flop. Their explicit forms in this case are 
shown in the appendix D. 

The densities ff h (pnp'n'; J£) and {^(pnp'n'; J n ) are supplied by the nuclear structure 
calculations. As an example, we discuss below the j3/3 decay 48 Ca — > A8 Ti. 

We first consider the case when the intermediate nucleus 48 Sc and the final nucleus 48 Ti 
are described, respectively, as one-particle one-hole and two-particle two-hole excitations 
on 48 Ca, i.e., 



\r a M) = E(H^> J ^ 

pn 



|0/> = E N(pp')N{nn')(p P 'nn';r\tf) [(aja£) r *(tw)i-j |0t), (49) 

with N{pp') = (1 + <5 P p/)~5. One gets from (|37|) 

(fipminn*, Jl) = (J:\p2n2) E E /(P3n3i^)(0;bpW / ;/ 7r }(-) ni+P3+J+/ 

p>p'n>n' pznsl'* 

x j ^ ^ j I A(PPOA(nn% pi 5 nni Vp 3 <W ( 50 ) 
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and 



p p c i(pin 1 p 2 n 2 ] J*) 



V I(Oj\pp'nn';r){-) ni+p i +J+I 



p>p'n>n' I* 



I>1 " [ J j. \ Pi(pp')Pi(nn')5 PPl 5 nni 5 p , P2 5 n , n2 



n 2 P2 



where 



Pj(pp f ) = N{pp') 



1 - 



\J+p+p' 



(p p') 



Within the CA one can use the closure relation 



J2fa n 3\ J a)( J a\P2 n 2) = <W 3 5 n 



2«3 ' 



(51) 



(52) 



(53) 



which leads from fl50D to fl5"l"|). 

On the other hand, within the QRPA formulation, and after solving the BCS equations 
for the intermediate nucleus 48 Sc f25j, the two-body density matrix becomes 



p ph (pnp'n';J^)= u n v p Xj,(pn) + u p v n Yj,(pn) u p >v n >Xj z (pn')+u nl v p >Yj^(pn) , (54) 



where all the notation has the standard meaning [|17], |25j • One should bear in mind that 
when the QRPA is used, the energies ujj* that appear in the matrix radial integrals are 
the solutions of the RPA problem and not the excitation energies of the intermediate 
nucleus relative to the initial nucleus. 



In particular, in the single mode model |fL6] , where there is only one intermediate 
state for each J 71 (and which seems to be a reasonable first order approximation for the 



(3(3 decays of 48 Ca and 100 Mo nuclei |EL E7|) 



p ph {pnpn; J") 



u p v n u n v p 



LO 



(55) 



where G{J n ) = G(pnpn; J n ), uj° = - [G(pppp; 0+) + G(nnnn; 0+)] /4, and = G(jj'jj'; J n ) 

are the particle-particle matrix elements. The intermediate states for Sc are: 

[O/7/2 (p) O/7/2 and the values of the ratios G{J + )/uiq for the 5 force can be found 



in table 1 of ref. [26| 
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5 Summarizing Discussion 

A straightforward derivation of the f3f3o u decay rate, based on the Fourier-Bessel ex- 
pansion of the transition amplitude, and the posterior application of the Racah algebra, 
has been performed without invoking the closure approximation. If necessary, this ap- 
proximation can be implemented, however, at any step of the calculation. It has been 
used for deriving the j3j3 0u formulas in refs. |], |], U, but not in refs. |L8|, |19| . 

To evaluate the nuclear matrix elements exhibited in eqs. (27) - (33) we only have to 
perform summations on the angular momenta and the intermediate virtual states. The 
successive terms rapidly decrease, because the radial integrals (f|5|) steadily diminish in 
magnitudes when the multipolarities L and L' are augmented |16|, [17]. The formulas 
become particularly simple when the harmonic oscillator basis is used. Then the Horie 
and Sasaki method |15| can be exploited for the evaluation of the radial form factors (f45[) , 



and the equations displayed in the appendix D can be used. 

The present formalism is especially suitable for the nuclear structure in which the 
summation on the intermediate states is unavoidable, such as the QRPA. The closure 
approximation then just connotes that the variation of the energy denominators with 
nuclear excitation is not considered. Evidently this does not lead to a major simplification 
in the numerical calculation. For example, because of p ph (pnp'n'; J£) given by (|54|), the 



summation in (|44| ) on different states a with the same J 71 " persists, although we do the 
replacement utjn — >< ujj* >= ujjtt. 

Contrarily, the use of the closure approximation is mandatory, and can be implemented 
easily as described in the last section, when the study is done in the shell-model framework, 
i.e., when one possesses information only on ff d , or equivalently on the + nuclear wave 
functions for the initial and final states. In this case the matrix element (HIT) reads 



M GT = -J2(-) L E P?(P"pV;J w )W' UJ (pn)W' L ij(pV)7ei £ (pnA> J *). 

LJ n pp'nn' 

The main difference between the formalism presented here and those published so far 
IH|, 0, |4|, 18, |19|] is its simplicity. As such it is more suitable for the numerical calculations. 



Let us underscore a few points in this regard: 

1) While in the neutrino potential formalisms |1], 0, |J one deals with two-body matrix 
elements, which lead to rather complicated analytic expressions for the (3(3q u moments, 
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we only have to handle the well known one-body operators (p6|). It could be illustrative 
to compare our result (|3T| ) for the matrix element Mr with eqs. (3.65) to (3.68) in the 
Tomoda's report H. 

2) At variance with the formalism developed by Vergados et al, ]TS[] , the results shown 
here are not limited to the employment of harmonic oscillator one-particle wave functions. 
Besides we totally avoid the usage of the Moshinsky-Brody transformation brackets, which 
now and then could be cumbersome. 

3) There are as well several substantial differences with the works of Suhonen et al. |n| , 
where the Fourier-Bessel expansion has also been used. First, they obtain different and 
more complex results for Mp and Mqt- Second, they do not exhibit the explicit structure 
of for the remaining matrix elements, given here by eqs. (29) - (33), but only show 
their general layout. Yet, this layout cannot be used for any practical purpose. Third, 
instead of dealing with the plain nuclear shell model, they operate in a relativistic quark 
confinement model. Fourth, their formulation is limited to the QRPA approximation as 
well as to the harmonic oscillator basis. 

In summary, we believe that the present formalism simplifies the nuclear structure 
evaluation of the (5f3o v matrix elements to a large extent. The formulation is applicable 
as well to matrix elements that appear in some supersymmetric contributions 
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Appendix A: Matrix elements Mx(k, n) 

After integrating on (ix and dy, as indicated in eq. the matrix elements (|11|) and 
(JlJ) read 

M_p(k, n) = gl( F \Y,T+^\N)(N\J2T+e-^\i), (A.l) 



M GT (k,7v) 



<?i(HE^>« eik ' r ">> • HE 7 "™^ ik ' r " l l 7 )' 

n rra 



(A.2) 



M^(k,iv) = -2^,(F|Vr+k-r n e jk - r l7v)(^|Vr+e- k - r -|/), (A.3) 



n 



m 



M' GT (k,N) = -2z^(F|^T+k.r^ n e ik - r ™|Ar).(iv|^T> me -*- r ™|7), (A.4) 



M K (k, jv) 



1 2M, 



A' 



I n 



r+<T n e ik - r "|iv)x 



(jv|^t+ p m e- lk - rm + e- 4k - rm p m + / w Vxcr m e- ik - rm |/) 



ikr,, 



(A.5) 



M T (k,N) = 5> n > n ® (k® r n )( 2 )]We ik - r "|iv) ■ (*| E?>^ k - r ™|/) 5 



v/3 



(A.6) 



M P (k,N) = -V2i^ v {V3<F|E^>„®(k®r^ 

(f| E r+cTne*"-^) • (iv| $>+(k® r m )( 1 )e- ik ' r - |i> } . (A.7) 

ra rra J 
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Appendix B: Electron matrix elements and phase-space 

factors 



The leptonic factors in eq. ( |T8| ) are 



Li{e 1 ,e 2 ) = (-1) 1/2 S '*x\>[g-M) - /i(ei)<7--p 1 ][/i(e 2 )<T-p 2 + ^_ 1 (e 2 )]x- s ^ 
L 2 (ei, e 2 ) = (ei - e 2 )(-l) 1/2 - fi 2 X t / j^ ^^^y^^^ . ^ + f^g^^ . Pl ] x _ s ,. 

^ 3 (ei,e 2 ) = ^(-l) 1/2+s ' 2 x; ; [^i(ei)^i(e 2 ) + / 1 (e 1 )/ 1 (e 2 )o--p l0 --p 2 ] X _ 4 (B.l) 

^(6 1; e 2 ) = ^(-l) 1/2+ ^4{[A(e 1 )/_ 1 (e 2 )+^ 1 ( ei )^i( e2 )]^-Pi 

- [/-i(ei)/i(e 2 ) + <7_ 1 (e 1 )<7i(e 2 )]<r • p 2 } X-s' 2 , 
L 6 {e u e 2 ) = ^(-l) 1/2+s ' 2 x; ; {[^i( e i)/-i( e2 ) + /-i(6 1 k- 1 (e 2 )] 

- [gi(ei) /i(e 2 ) + A(ei)#i(e 2 )]cr • Pl cr ■ p 2 } 

where all the notation has the usual meaning ||. 

The electron phase-space factors J-^{Tq) that appear in eq. fl2"4|) are: 

^i(To) = T (30 + 60T + 40T 2 + 10T 3 + T 4 )/30, 

^ 2 (T ) = T 4 (70 + 77T + 14T 2 + T 3 )/420, 

F 3 (T ) = T 3 (10 + 10T + T 2 )/30, 

J- 4 (T ) = r Q 2 (30 + 35T + 10r o 2 + T Q 3 )/135, 

^5(^0) = r [60To + 80r 2 + 30r 3 + 3r o 4 + e(60 + 90To + 40T 2 + 5T 3 )]/45, (B.2) 

J- 6 (T ) = 2To(12 + 18T + 8T 2 + T 3 )/(3R), 

F 7 (T ) = 4T [60T + 100T 2 + 55T 3 + 12T 4 + T 5 + £(60 + 90T + 45T 2 + 10T 3 + T 4 )]/(45R), 

JF 8 (T Q ) = T [100T 2 + 150T 3 + 73T 4 + 14T 5 + T 6 + 2£(60T + 100T 2 + 55T 3 + 12T 4 + T 5 ) 

+ £ 2 (60 + 90T + 45T 2 + 10T 3 + T 4 )]/135, 

J- 9 (T ) = 4To(60 + 90To + 45T 2 + 10T 3 + T 4 )/(15R 2 ), 



with 
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Appendix C: Derivation of the final formulas for the nuclear 

moments 



Below we give the details on the derivation of the last term in (|31|). First we rewrite 
flOQ as 

M*(k,*) = -*^k.(F|£r> n e* k - r >>x 
ZM N n 

M E ^e- ik ' rm [2 Pm - k + if w <r m xk] \i), (C.l) 



and then we express the vector product, involving the nucleon momentum term in ( |C.1| ), 
in spherical coordinates 



(*1 J2r^e ik ^\ N )(N\ J2T+e-^p-»'\i). (C.2) 



rn 



After performing the multipole expansion fl2"5|) and handling some straightforward Racah 
algebra, we obtain 

M<*\k, N )=*^VE(A*)*k £ f M(L010|.0) 

1VlN LL'JJ'M'MW 

(i, -vk P \j', -m')(i, -z/L'MVM)r Kp (k)r;, M ,(k)(F|s LLJ , M ,(fc)|iv)(iv|p LV M 

(C.3) 

where the tensor operators S° LLJ ,(k) and Pz/j(&) are defined in (|26|) . Finally, the angular 
integration allows us to perform the summations on the angular momentum projections 
and obtain 

/^ k M(P)(k,iv) = 5j^v/6(47r) 2 Jfc E i L ~ L 'L | \ \ M(L010|L'0) 

WS^j(*)k>-HPi/j(*)W- (C4) 
This result, together with the eq. (|i~4"|) , yields the last term in eq. ([31]). 
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Appendix D: Radial form factors for the harmonic oscillator 

wave functions 



Following the Horie and Sasaki method (15| the radial integral ( fj5|) can be expressed 



as: 



T^LL'(P n P n ''^Js) = [M(p,n)M(p',n')] 1/2 ^ a m (p,n)a m <(p f ,n)fZ L ,(m,m'-,ujjn), 

mm' 

(D.l) 

where 

M{n p l p ,n n l n ) = 2^ +n -n p \n n \(2l p + 2n p + l)\\(2l n + 2n n + l)\\, (D.2) 



n p \l n n \ (2l p + 2n p + l)\\(2l n + 2n n + l)\\ 
k £f =s \ k J { k' ) (2l p + 2k + l)\\(2l n + 2k' + l)\\ • 



a>i p +i n +2s(n p l p ,nj n ) = x 1 "'" 11 "" 1 — 



(D.3) 



and 



f£ L ,(m,m';u JS ) = J> 2 J ^—^L, ^— -L' ) J£(ujj s ), (DA) 



J£(wj S ) = (2v)-' t Rj dkk w+K e~ k / 2u v{k;ujj s ). (D.5) 
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